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Weak coupling of internal degrees of freedom of molecules to the translational degrees of freedom 
of a fluid results in additional modes of motion for density fluctuations. These new modes affect the 
spectral distribution of light scattered by density fluctuations so that the Landau-Placzek ratio is not 
satisfied. The case of thermal relaxation with a single relaxation time is worked out in detail. For- 
mulas for the spectral distribution of the scattered light, for the ratio of the intensities of the central 
(Rayleigh) to the Brillouin components and for the phonon velocity are derived and applied to carbon 
disulfide and carbon tetrachloride. The data for carbon tetrachloride are shown to be inconsistent 
with the single relaxation time model for thermal relaxation. 
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1. Introduction 

In this paper we are concerned with density fluctuations in a fluid in which internal degrees 
of freedom of the molecules are weakly coupled to the translational degrees of freedom of the fluid. 
Thermal relaxation is an example of the type of processes we have in mind. Problems of this 
type are of interest because it is now possible to experimentally probe the frequency spectrum of 
density fluctuations with light scattering experiments using a laser as the light source [l]. 1 
The coupling of internal degrees of freedom to the translational motion means that the decay of 
density fluctuations will proceed by more ways than the three "normal modes" usually considered. 
This is reflected in the spectral distribution of light scattered by density fluctuations. In this 
paper we investigate a relatively simple case involving thermal relaxation with a single relaxation 
time. The existence of the additional mode is shown and the effect of this mode on the spectrum 
of the scattered light is examined in detail. 

Light scattering experiments provide a Fourier analysis of the density variations in a fluid. 
The variation of the intensity of the scattered light with the scattering angle 6 (and therefore with 
the change in the wave vector of the scattered light) Fourier decompose the spatial dependence 
of the fluctuations while the shifts in the frequency of the scattered light decompose the time 
dependence of the fluctuations [2]. Such experiments enable us to study collective motions in 
the fluid without seriously disturbing the fluid. On the other hand, it is also possible to use light 
scattering measurements in conjunction with a model for the fluctuations to measure several of 
the bulk properties of the fluid. For example, the sound velocity and the sound absorption co- 
efficient can be obtained by measuring the shift in frequency and the width of the Brillouin lines. 

In some circumstances it is possible to obtain from light scattering experiments information 
about the structure of the correlation functions whose time integrals are the transport coefficients 
[3]. An example of this type provides the basis for the calculations we present. 

This paper consists of three parts. In the first part we review the relationship between density 
fluctuations and light scattering. The phenomenological approach of Einstein and Smoluchow- 



1 Figures in brackets indicate the literature references at the end of this paper. 
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ski [4] is discussed and compared with the formulation of scattering problems of Van Hove [5], 
The generalized structure factor S{k, o>) is introduced and some of its properties are described. 

The second part of the paper consists of a detailed analysis of the frequency spectrum of 
density fluctuations in a fluid in which thermal relaxation of vibrational degrees of freedom can 
occur. Particular emphasis is placed on single relaxation time processes. The results for 
single relaxation time processes are summarized at the end so that readers not wishing to go 
through the details of the calculation will not lose heart. 

The third part of the paper consists of a comparison of available experimental data with 
the analysis of the density fluctuations. The primary example is concerned with CS2, a substance 
whose internal degrees of freedom relax with a single relaxation time. We also demonstrate 
that the data for CCI4 are not in agreement with the predictions for a single relaxation time. 



2. Light Scattering Formalism 
2.1. Phenomenological Theory of Einstein-Smoluchowski 

The random thermal motion of molecules in a, fluid produces fluctuations in the density and 
also in the orientation of molecules in volumes small compared to the wavelength of the incident 
light. These fluctuations result in local variations in the dielectric constant and therefore scatter 
light. We are concerned in this paper only with fluctuations in the density. Orientation fluc- 
tuations result in the depolarization of the scattered light making it possible to experimentally 
separate the scattering by density fluctuations, which is fully polarized, from the scattering by 
orientation fluctuations [6]. 

The intensity of the scattered light is 



/(R,o>) = /o 



Nkf 
16tt 2 R 



Jsin 2 ^([€(A % co)] 2 ). (1) 



In eq (1) incident plane polarized light of intensity h and wave vector k* is scattered at the origin 
and is observed at R. There are N molecules in the scattering volume. The angle between the 
electric vector of the incident wave and R is <p; e(k, co) is the Fourier component of the fluctuation 
in the dielectric constant. The shift in the angular frequency of the scattered light is a) and the 
change in the wave vector in the medium of the scattered light is k; since only the direction of 
k changes 

k = 2nh sin 0/2. (2) 

The index of refraction of the scattering fluid is n and the scattering angle is 6. The angular 
brackets (...) indicate an ensemble average over the initial states of the system. 

Direct calculation of e(k, a>) is avoided by assuming that fluctuations in the dielectric constant 
are due to fluctuations in the density and the temperature; 

Ae = (deldp) T Ap + (de/dT) p kT. (3) 

The contribution of the temperature fluctuations is ignored; we assume that (de/dp)r ^(dejdT)p. 
Equation (1) is now reduced to 



7(R,o>) = /o 



Nkf 



l6ir 2 R 2 



sin 2 <p(deldp) 2 T ([p(k,co)Y) (4) 



where p(/r, a>) is a Fourier component of the density fluctuation. The problem is now one of 
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calculating p(A, co) from the appropriate kinetic model of the fluid. For dense fluids the appro- 
priate model is described by the linearized hydrodynamic equations of irreversible thermodynamics. 
We should note that (de/dp) T has been assumed to be independent of the shift in the frequency 
of the scattered light, co. If this is not the case it is probably also true that the separation of 
Ae into a thermodynamic derivative times a fluctuation term is not a meaningful procedure. Of 
course the value of (de/dp)r may vary somewhat as the frequency of the incident light changes, 
reflecting the frequency dependence of the dielectric constant. 



2.2. Molecular Theory 

Before we undertake the calculation of p(k\ co) we review the light scattering formalism ob- 
tained by adapting Van Hove's neutron scattering paper [5] to light scattering. Komarov and 
Fisher [7] have shown that the intensity of light scattered from a fluid of N molecules of effective 
polarizability a is 



/(R, o>) = /o 



Na 2 kj 



2ttR 2 



sin 2 <pS(k\ 0)). (5) 



The important change from eq (4) is to replace the mean square fluctuation ([p{k\ co)] 2 ) by the 
generalized structure factor S(k\ co) which is the space and time Fourier transform of the two-body 
correlation function. The correlation function is defined by Van Hove to be 

C(r, t) = N~* /J f dr'8[r + rKO)-p']8[p'-pj(0]y (6) 

For long times and sufficiently large r, G(i\ /) reduces to the autocorrelated density 

G'(r, t) = N~ 1 j c/r'<p[r'-r(0), 0]p(r\ *)>■ (7) 

Equation (7) is appropriate to light scattering in fluids [8|. Care must be taken to use 

S(k\ oj)= dr I dte ik r e - i(ot G(r, t) (8) 

only to describe the fully polarized part of the scattered light. The inclusion of angular correla- 
tions, which result in depolarization, is a more complicated problem than the one we consider 
here [9|. 

In this paper we are concerned with §{k :, co) as defined by eqs (7) and (8); 

S(k H (o)=(p(k-, co)p(-A)). (9) 

Thus eqs (4) and (5) predict the same frequency distribution for light scattered by density fluctua- 
tions. A useful property of S(k\ co) is the sum rule 

S(k) = j- \~ dw£(k, co)=<p(A)p(-A)>. (10) 

S{k) is the ordinary structure factor. Finally, we note that c£(A , co) is an even function of a) at high 
temperatures; that is Hoj/k'nT <^ 1 where kn is Boltzman's constant and T is the absolute tempera- 
ture. In this paper we assume that the inequality is satisfied. 
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3. Density Fluctuations When the Internal Modes Decay With a Single Relaxation 

Time 

We are concerned with the calculation of £(k, co) for a fluid whose molecules possess internal 
degrees of freedom which are weakly coupled to the translational modes of the fluid. Further 
we shall suppose that the transfer of energy from the internal degrees to the translational degrees 
of freedom is described by a single relaxation time process [10, 11]. 

The procedure to be used is to solve the linearized hydrodynamic equations for p(k, co) in 
terms of an initial fluctuation p(k) [5]. The use of initial values facilitates calculation of the aver- 
age over initial states indicated by (. . .). 

An equivalent procedure is the hydrodynamic fluctuation theory of Landau and Lifshitz [12]. 
This method has been used by Rytov to discuss fluctuations in a viscoelastic medium [13]. 

The linearized hydrodynamic equations for the system are the continuity equation 

dpi/dt + po div v = 0, (11) 

the longitudinal part of the Navier-Stokes equation (suitably modified to allow for a frequency 
dependent bulk viscosity) 

+ (f * + 1 fc) S rad diV ▼+[%'(*-*') S rad div V iP'Wn (12) 

and the energy transport equation. 

PoCvidTJdt) - [c v (y-l)/p] (dpildt) - XV 2 7\ = 0. (13) 

Here p = p + pi is the number density, T=T {) +Ti is the temperature; p and T being the equi- 
librium values. The shear viscosity is rjs, the bulk or volume viscosity consists of a frequency 
independent part i) v and frequency dependent part which is the Fourier transform of t)\t). The 
low frequency (adiabatic) sound speed is c , the thermal expansion coefficient is j8, the thermal 
conductivity is X and the ratio of the specific heat at constant pressure c p to the specific heat at 
constant volume c v is denoted by y. 

The analysis proceeds as in reference 8. First the Fourier (space) and Laplace (time) trans- 
forms are taken of eqs (11), (12), and (13). Then we solve for 



p(&, s)= dire-*-* dte~ 8t pi(r, t) 

Jv J0 



(14) 



in terms of the initial value 



The result is 



p(k)= I dre-*r p ( r , 0) (15) 

Jv 



P(A) G(s) (l0) 
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and 



F( S ) = s 2 + s [ak 2 + b k 2 + b'(s)k 2 ] + ab k 4 + ab'(s)k A + c 2 k 2 (l - 1/y) 



G(s) = s 3 + s 2 [ak 2 + bok 2 + b'(s)k 2 ] + s [c 2 k 2 + ab k A + ab'(s)k A ] + c 2 ak 4 ly. 



(17) 



(18) 



The quantities 



a = k/p {) c v 

b = (4>l3r) 8 + r) n )/p Q 

1 f x 
b'(s) = — e- st r]\t)dt 
po Jo 



(19) 



have been introduced to simplify t lie notation. Equation 19 for b'(s) applies to any frequency 
dependent bulk viscosity. We will be concerned with b'(s) for a single relaxation time process; 



b'(s) = b,l(l+ST) 



(20) 



where r is the relaxation time. 
The dispersion equation is 

G(s) = (21) 

and is used in ultrasonics with s replaced by ico. It is known from ultrasonics [10] that 

6i = (c2-c?)t (22) 

where c x is the infinite frequency sound speed. In thermal relaxation [11) 



\Cp — Cy)Cl 

(c v — ci)c p 



Col", 



where c/ is the specific heat of the vibrational degrees of freedom. 
To compute S{k, a>) we first observe that 



Jo 



p(k, co) = 2 Re dte-^pik, t) 
Jo 

where p(A\ t) is the inverse Laplace transform of p(k, s). It follows that 

£(k, (o) = (p(k)p(-k)) {2 Re F(ioj)/G(io))} 
= (p(k)p(-k))a(k,a>). 



(23) 



(24) 



An exact expression for cr(/r, o>) may be readily obtained by replacing s by ia) in eq (16) and taking 
twice the real part. The resulting expression is quite involved; direct substitution leads to 



a(k, <o) = 2[N 1 D 1 + N 2 D 2 ]l[Dl + I%] 



(25) 
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Figure 1. The central component of cr(k, oj) for scattering FIGURE 2. The phonon component of cr(k, a>) for scattering 
in CS 2 at 20 °C as a function of oj. in CS 2 at 20 °C as a function of oj. 

When a He-Ne laser is the source, A=8.608X 1() 4 cm"' corresponds to seat- When a He-Ne laser is the source A =8.608 X 1() 4 cm" 1 corresponds to scat- 

tering at 30° terin S at 30° 



where 



and 



Ni =- to 2 + ab k 4 + c 2 k 2 (l - 1/y) + {ab,k 4 + b x k 2 oj 2 r)l{\ + coV), 

N 2 = co[ak 2 -h 6 A- 2 + {b x k 2 - ab x k*7)l{\ + co 2 r 2 )] , 

£>! =- o) 2 (aA' 2 + W 2 ) + ciak^ly + iab^c^T-o^b^lil + w 2 t 2 ) 

D 2 = a>[-a) 2 + C& 2 + abok 4 + (/>!yc 2 w 2 r + a6^ 4 )/(l + o> 2 r 2 )] . 



(26) 



(27) 



Now we have an exact expression for the frequency distribution of the scattered light. If 
we wish to be able to interpret this in terms of the properties of the material, we must develop 
a way to pick out the significant portions of cr(A, co) for different values of co. We have evaluated 
cr(A, a)) for a representative set of material parameters and A = 8.608 X 10 4 ; the results are shown 
in figure 1 and figure 2. This corresponds to = 30° for scattering in CS 2 using a He-Ne laser. 
To obtain as much useful information as possible from such a pattern it is necessary to write eq 
(25) as a sum of terms which are individually important only over a restricted range of frequencies 
co. The denominator does not obviously factor so an indirect approach is needed. One method 
is to pick out the important terms, say for small a> and to discard the rest. Another approach, 
one used profitably in reference 8, is to approximately compute the inverse Laplace transform 
of p{k, s) and then to compute p(A:, o>) by means of eq (23). The virtue of this method is that alge- 
braic expressions for the modes of motion of the density fluctuations are obtained. The difficulty 
of course is to obtain good approximate solutions to the dispersion equation, eq (21). We shall 
make use of both approaches to investigate the properties of cr(/r, a>). 



4. Properties of cr(A\ co) 

In this section we are concerned with the construction of a good approximation to a(k, co) 
which will permit interpretation of spectral distribution curves in terms of the properties of the 
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scattering fluid. The results of this section are contained in eqs (43), (45), and (46), and are listed, 
for convenience, in subsection 4.6. 



4.1. Approximate Solution of the Dispersion Equation 

When b'(s) is given by eq (20), the transform for the density, eq (16), is 

p(k,s) = Fi(s) 
p(k) G t (s) 

where 

Fi(s) = (1 + st)s 2 + s[(1 + sr)(ak 2 + b {) k 2 ) + b,k 2 ] + (1 + sr) [ab k 4 + cjA^l - 1/y)] + abxk 4 (28) 

and 

ftto = rs 4 + s 3 [1 + r(ak 2 + bk 2 )] + s 2 [ak 2 + b Q k 2 + b x k 2 + r(c?A: 2 + a/U 4 )] 

+ s [elk 2 + aio* 4 + abtk* + r(c 2 ak A /y)] + c 2 ak 4 ly. (29) 

The first step in computing p(k, t) is to obtain approximate roots of 

Gi(s) = 0. (30) 

It is convenient to introduce the dimensionless quantities 

Y= s/cok 

a = ak 2 /cok 

po = b k 2 lcok (31) 

ft = bttf/cok 

D = cokr 

into eq (30) with the result 

DY 4 + F 3 [l + D (a + j3<>)] + Y\a ^p o +Pi^D(l + ap )]^Y[l^apo+apv-\-Daly]-\-aly = 0. 

(32) 

For interesting values of the change in the wave vector, k, 

D~l 9 /8i«l 
while 
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First we look for solutions of eq (31) which are of order unity. The relevant parts of eq (31) are 

DY" + Y s + (D + ^)Y 2 + Y^0; (33) 

we have neglected terms of order a and a 2 . One root of eq (33) is zero so we also must examine 
eq (32) for solutions of order a. The dominant terms are then 

Y+a/y^O. (34) 

Here terms of order a 2 and smaller have been dropped. The solution to eq (34) is 

Y=-aly (35) 

or 

s=—\k 2 lp(,Cp 

which also appears in classical theory. 

When the F=0 term is removed from eq (33) we have 



Y S + Y 2 ID + 



D + /3, 
D 



F + i = 0. 



(36) 



The formal factoring of eq (36) is 

(Y+A)(Y 2 + BY+C) = 0. (37) 

Two solutions are 

Y=-BI2±iC l i 2 [l-B 2 l4C] 1 ' 2 (38) 

which correspond to the phonon modes. The third solution, 

Y=-A 

yields a second nonpropagating mode. In order that damping occur, it is necessary that A > 
be satisfied. For phonons to exist it is necessary that Z? 2 /4C ^ 1 so 



Y^-BI2±iO' 2 . 
Comparison of eqs (36) and (37) shows that 

AC=l/D 
therefore 

Y=-HCD. 
Now multiply eqs (38') and (39') by cok to obtain 

s = T B ± ivk; vk — CokC 112 
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(38') 



(39') 



(40) 



and 



S — — (%\v 2 T. 

The lifetime of the phonons is (T^) -1 and their speed is v. These are as yet unspecified quantities 
to be extracted from the exact expression of cr(A , a>). Before we do this we shall obtain an approxi- 
mate formula for <j{k :, to). 

4.2. Calculation of cr(A\ co) 

We take &(s) to be 

Gi(s) = r(s + ak 2 /y) (s + T B ~ ivk) (s + T B + ivk) (s + 4/v 2 t). (41) 

Combining eqs (28) and (41) it is a straight forward process to obtain the inverse Laplace trans- 
form of p(k, 5). Ignoring small terms, we obtain 



P(*, t) 
p(k) ' 



= (l-l/y)e-^7poC,+ 



+ 



(4 - cg)A* - (v 2 lcj - 1) (4Iv*t 2 + c§A*(l -1/y)) " 
cj/tfV + ^ 2 A 2 



e — cUlv 2 r 



[1-cfMl - 1/y)] (v*k* + <%IW)-(<*,--<%W 
z; 2 A 2 + cJ/vV 



X e-'V cos (vkt). (42) 



Next we apply the operation indicated in eqs (23) and (24) to obtain 

2AA- 2 /poCy, 



<r\k % co) -(1- 1/y) 



(XP/poC p ) 2 + co 2 



+ 



(c 2 x - c 2 )A 2 - (v 2 lc 2 - \)(ctlv*T 2 + c 2 A- 2 (l - 1/y)) 
c?>/z; 4 t 2 + ?; 2 P 



2cllv 2 r 



+ 



c^V-fco 2 
[1 -c 2 Ml - 1/y)] [i; 2 A: 2 + c 2 /z; 2 t 2 ] -(cj-c 2 )k 2 



CoJv 4 t 2 + 1; 2 A -2 



+ : 



r« 



iri+(a)-vk) 2 r 2 + ( W + 



»*)»J" 



(43) 



The prime means that <x' is an approximation to (r(k, co), which is given by eq (25). 

The first term in eq (43) corresponds to decay of a fluctuation by a thermal diffusion process. 
We refer to it as the thermal mode. The second term also represents a nonpropagating type of 
decay which is coupled to the internal degrees of freedom of the molecules. The last term rep- 
resents the phonon modes. Equation (43) is derived on the assumption that (aA 2 /y)r is much 
less than one. If this is not the case, eq (43) is a poor representation of the nonpropagating modes. 



4.3. Determination of v{k) and T B 

Thus far the phonon speed v and width T B are undetermined quantities. In this subsection 
we remedy this by comparing the phonon terms of eq (43) with the significant parts of or(k, co), 
eq (25), when co ~ vk. 
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Figure 3. The phonon velocity v(k) as a function of the 
wave vector kfor CS 2 at 20° C. 

The experimental point (x) is taken from reference 18. Equation (45) and 
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The direct evaluation of a(k, co) shows that the phonon peak coincides with the vanishing of 
the imaginary part of the dispersion relation, 



A>(co) = 0. 



(44) 



The positive root of eq (44) determines the phonon frequency v{k)k. Neglecting the term involving 
the thermal diffusivity, a, eq (44) is 



- co 2 + cgA? + btftfrUX + <o 2 t 2 ) = 0. 



The solution is 



(cor) 2 = (v(k)kr) 2 = 1/2[(cooAt) 2 - 1] + 1/2[(1 -c 2 c ^ 2 t 2 ) 2 + 4c 2 ^ 2 t 2 ] 1 / 2 . 



(45) 



From eq (45) we see that as k— » 0, v(k)—>'Co and as k— *oo, v(k) — > c x . The variation of v(k) with 
k is indicated in figure 3 

The width T B is obtained by observing that if D 2 (v(k)k) = 0, D 1 (v(k)k) = -2[v(k)kfT B . The 
minus sign is necessary so that Y B > be satisfied. Direct substitution in eq (26) yields 



2T B = ak 2 + b k 2 - 



c|\ (aP\ 

v 2 A 7 ) 



b,k 2 



l+V 2 k 2 T 2 



(l-ak 2 r) 



(46) 



It should be noted that the width due to the relaxation is not simply added to the classical absorp- 
tion term ak 2 (\ — 1/y) -\- b k 2 . In practice this difference may not be significant, although for 
CS2 it amounts to about 5 percent of T B . 

The method used to obtain v(k) and T B should be applicable to more complicated situations 
although it would be advisable to examine Di((o) and D 2 (co) numerically and to verify that D 2 (v(k)k)=0 
for each new situation. 

4.4. Comparison With the Exact Expression 

The approximate frequency distribution contained in eq (43) has been compared numerically 
with the exact expression, eq (25). This is illustrated in figure 4 where the percent deviation of the 
approximate expression from the exact one is shown as a function of frequency. The parameters 
are the ones used to obtain figure 1. The deviation in the central components is probably due to 
a small error in the width of the thermal diffusion mode. A decrease of 1 percent in that width 
would eliminate most of the difference between cr{k, co) and cr'(k, co). The deviation in the phonon 
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component indicates that the maxima of the exact and approximate expressions do not quite coin- 
cide. This results in the oscillating deviation shown in figure 4. Deviations greater than 1 per- 
cent occur when the magnitude of cr'{k, co) has fallen to less than 1 percent of its maximum value 
o"'(/r, 0). The sum rule, eq (10), is satisfied by the approximate expression cr'(k, co) to within 1 
percent using this set of parameters. 

It is, of course, possible to improve the accuracy of <x'(A, co) by obtaining more accurate solu- 
tions to the dispersion equation. This would involve using the solutions we have found, eqs (35) 
and (40), as the starting point of an iteration of the dispersion equation. The resulting formula 
for o~'(A% co) would be much more complicated than eq (43). 

4.5. Intensity Ratio 

The ratio of the intensity of the unshifted (central) components of the scattered light to the 
intensity of the Brillouin components is a quantity which is readily obtained experimentally. 
This ratio, J> =I c I2Ib, is easily obtained from eq (43) by integrating the individual terms: 



>=- 



(cj - eg)* 2 ~ (v 2 lcl - 1) [(4Iv*t 2 ) + clk\l - 1/y)] 

C 4 IV 4 T 2 + V 2 k 2 

[1 - dlvM - 1/y)] [v 2 k 2 + c 2 Jv 2 r 2 ] - {cl - c 2 )k 2 



(l-l/y)+ 



(47) 



C 4 IV 4 T 2 + V 2 k 2 

At low phonon frequencies vkr < 1 this reduces to the Landau-Placzek result 

J? = y-1. (48) 

At large phonon frequencies (vkr > 1) a more involved expression than eq (48) is obtained; 

^ = (y-l)[\ + {yl(y-\)}{(cl-cl)lcl}]. (49) 

The high frequency limit, eq (49) was obtained earlier by Rytov [13]. 
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The ratio of intensities J^, given by eq (47) agrees within 1 percent of J> obtained by numerical 
integration of the exact expression for cr(k, co). From an experimental point of view, eq (47) should 
be taken to be an upper bound on J> . This is because it may be difficult to detect all of the light 
scattered by the second nonpropagating mode; the difficulty arises from the large width of this 
component. 



4.6. Summary of the Properties of cr(k, a>) 

In the previous paragraphs we have shown that cr(/r, co), the frequency distribution function, 
may be approximately represented as a sum of four Lorentzian curves; 



cr'(*, ©) = (1- 1/y): 



2\k 2 lp Q c l , 



(kk*lpoc p r + a> 2 



+ 



(c 2 ~ cl) k 2 - (v*ld - 1) (ctlv 4 T* + c 2 k 2 (l - 1/y) ) 



ctlv A r 2 + v 2 k 2 



2c 2 lv 2 r 



(cllv 2 T) 2 "f 0) 2 _ 



(43) 



[I-c 2 ) lv 2 (l-lly)][v i k 2 + c 2 lv 2 T 2 ]-(cl-c 2 )k 2 



ctlv 4 T 2 + V 2 k 2 



Tb 



: + 



r B 



Tl + (C0 - vk) 2 ri + (co + vk) 2 _ 



(44) 



where the phonon speed v(k) satisfies 

[v(k)kT] 2 = 1/2 [ (Coc krf - 1] + 1/2 [ (1 - c 2 J 2 r 2 ) 2 + 4c 2 £- 2 t 2 ] l > 2 



(45) 



and the phonon width is 



2r«=aF + M 2 -@V— 1 



l+[v(k)kT)] 



(I-oFt). 



(46) 



The ratio of the central components to that of the Brillouin lines is, in the high frequency limit 
(v(k)-> c x ) 



if = (y - 1)[1 + {y/(y - 1)} {(cl - c 2 )/c 2 }]. 



(49) 



5. Comparison With Experiment 
5.1. CS 2 

In the previous sections, various points have been illustrated by using parameters appropriate 
to carbon disulfide (CS2). For example, Figures 1, 2, and 3 show a (/c, co) and v (k) for CS 2 as 
predicted by egs (25) and (45). Comparison with experiment is limited to the phonon speed and 
to t/, the ratio of the intensities of the central and Brillouin components. 

The parameters for CS 2 at approximately 20 °C are listed in table 1. Several measurements 
of the phonon speed are listed in table 2. To within the precision with which r and k are known, 
eq (45) predicts that v (k) lies within the experimental uncertainty in each case. 
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Table 1. Material parameters: CS-> at 20 °C 



TABLE 2. Phonon velocity and intensity ratio; CS2 



Parameter 


Value 


Rel 


erence 


a 


1.95 x 10- :) cm 2 /sec 




6 


b 


3.84 x 10- ;i cm 2 /sec 




14 


&i 


4.76 cm 2 /sec 




10 


Co 


1.142 x 10' cm/sec 




10 


T 


2.04 X 10 ! » sec 




10 


n 


1.675 




6 


p<> 


1.262 gm/cm 3 




6 


7 


1 .55 




10 





Experimental 


Calculated 










Velocity 


./ 


Velocity 


J 




m/ sec 




mlsec 




15 


1235 


— 


1241 


0. 83 


16 


1241 


■0.77 


1241 


.83 


17 


1250 


.71 


1241 


.83 


6 


1265 


. 93 


1241 


.83 


18 


1223 


— 


1212 


.69 



The experimental velocities have been reduced to 20 C using dv/dT—~S.2 m/sec- 
deg |6|. The velocities were calculated using eq (45). J willies were calculated using 
eq (47). The material parameters were taken from table I and k, the wave vector, was 
determined from information in the references. 



Three of the references contain measurements of J> . Here the agreement is not as close. 
The value of,/ predicted by eq (49) is 0.83. That the agreement with the phonon velocity is better 
than the agreement with J^ is not surprising. The limited sensitivity of the detectors used in these 
measurements makes it difficult to obtain accurate determinations of J> . The presence of the 
"second" nonpropagating mode aggravates the situation; unless the detector is quite sensitive, 
much of the light in that mode would be lost in the noise of the detector. 

5.2. CCI 4 

Sufficient data exist for carbon tetrachloride (CCI4) to warrant comparison with the prediction 
of the single relaxation time theory. This comparison indicates that only a part of the energy of 
the internal degrees of freedom relaxes with a time on the order of 10~ 10 — 10 n sec. Presumably 
the remainder of the energy relaxes more rapidly than can be detected by light scattering 
experiments. 

The parameters for CC1 4 at approximately 20 °C are listed in table 3. Four measurements 
of the phonon speed obtained from Brillouin scattering measurements and one acoustic measure- 
ment of the speed of sound in CCU are listed in table 4. The calculated velocities in table 4 are 
based on the assumption that only 75 percent of the internal specific heat is involved in the relaxa- 
tion process and that 7 = 5.12 X 10~ n sec. This value for r is consistent with the acoustic measure- 
ment. The agreement between the measured and calculated values of J^ is poor. This is not 
unexpected if another relaxation time exists for CCI4. This would imply the existence of another 
nonpropagating mode which would make J> larger than the prediction of eq (47). 

Also, for CCI4, there is a measurement of the width of the Brillouin line. Experimentally 
r w ~ 4 X 10 9 rad/sec while eq (46) predicts Y n ~ 5 X 10 9 rad/sec if 3/4 of the internal specific heat 
is involved in the relaxation process. 



TABLE 4. Phonon velocity and intensity ratio; CC1 4 



TABLE 3. Material parameters: CC1 4 at 20 °C 



Parameter 


Value 


Reference 


a 


1.78 x K) 3 cm 2 /sec 


6 


60 


5.98 X 10 •' enr/sec 


14 


6, 


3.59 x 10- 2 em 2 /sec 


10 


<"(! 


9.20 x 10 4 cm/sec 


10 


T 


5.12 X 10 " sec 




// 


1.472 


6 


P<> 


1.594 gm/cm 1 


6 


y 


1.46 


10 



The relaxation time t was estimated from acoustic measurements [19]. 





Experimental 


Calculated 










Velocity 


./ 


Velocity 


J* 




mfsec 




m/ sec 




6 


1040 


1. 10 


1024 


0.77 


16 


988 


0.75 


990 


.62 


17 


1000 


.72 


990 


.62 


20 


1002 


— 


990 


.62 


19 


972- 


acoustic- 




measurement 



The experimental velocities have been reduced to 20 °C using dvldT= — 3A m/sec — deg 
[6]. The velocities were calculated using eq (45). J? values were calculated using eq 
(47). The material parameters were taken from table 3 and k, the wave vector, was 
determined from information in the references. 
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In a previous publication [21] we stated that more than one relaxation time exists in CC1 4 
and therefore the single relaxation time theory did not apply. Although the available data are 
insufficient to fully verify this statement, the phonon speed comparison suggests that this state- 
ment is correct. 
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